Dynamics of photoexcited states in one-dimensional dimerized Mott insulators 
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Dynamical properties of photoexcited states are theoretically studied in a one-dimensional Mott 
insulator dimerized by the spin-Peierls instability. Numerical calculations combined with a per- 
turbative analysis have revealed that the lowest photoexcited state without nearest-neighbor in- 
teraction corresponds to an interdimer charge transfer excitation that belongs to dispersive exci- 
tations. This excited state destabilizes the dimerized phase, leading to a photoinduced inverse 
spin-Peierls transition. We discuss the purely electronic origin of midgap states that are observed 
in a latest photoexcitation experiment of an organic spin-Peierls compound, K-TCNQ (potassium- 
tetracyanoquinodimethane) . 

PACS numbers: 78.20.Bh, 71.35.-y, 71.10.Fd 



I. INTRODUCTION 

Photoirradiation is one of important methods to 
control physical properties of strongly correlated elec- 
tron systems iSi^i^^SiLSiSiiSiii Drastic changes of macro- 
scopic phases triggered by the photoirradiation, often 
called "photoinduced phase transitions (PIPTs)" ; 12 i 13 
are observed in a variety of materials where electrons 
are coupled with lattice systems A^LSiSiiSi 

Many PIPTs are understood in the framework of the 
adiabatic approximation. The key quantities are adia- 
batic potentials of few electronic states as a function of 
a lattice configuration (represented by Q). The adia- 
batic potential of the photoexcited electronic state has 
an energy minimum at Q — Q^ lin that is generally very 
different from the stable lattice configuration Q 9 min in the 
adiabatic potential of the ground state. According to the 
Franck-Condon theorem, the lattice configuration of the 
photoexcited state is initially at Q g min ( i.e., immediately 
after the photoexcitation from the ground state). Then, 
the lattice configuration moves from Q 9 min to Qf nin along 
the adiabatic potential of the photoexcited state. Such 
a photoinduced transformation of the lattice configura- 
tion takes place macroscopically in a PIPT for electron 
systems coupled with lattice systems ,I2iiii£ii£i 

The understanding of why the photoexcited state has 
a stable point different from that in the ground state is 
essential to reveal the fundamental physics of the PIPTs. 
To address this problem theoretically, we should first 
clarify how the photoexcitation induces charge transfer 
(CT) processes, although it is generally difficult to carry 
out in strongly correlated electron systemsiiiiSii2iS& 

As a typical example of strongly correlated electron 
systems coupled with lattice systems, we focus on a one- 
dimensional (ID) electron system with a dimerized lat- 
tice distortion, which is realized in an organic material, 
potassium-tetracyanoquinodimcthane (K-TCNQ). This 
material is a ID dimerized Mott insulator below the spin- 
Peierls transition temperature, T s p = 395Ki2i*22i2 3 i In 
1991, Koshihara et al. have demonstrated that an irradi- 



ation of a pulsed laser weakens the lattice dimerization of 
K-TCNQ in the dimerized phase, which is regarded as a 
photoinduced inverse spin-Pcicrls transition. 10 More re- 
cently, a high-resolution experiment has revealed several 
ultra-fast phenomena immediately after the photoirradi- 
ation, for example, the emergence of a midgap state in 
the reflectivity spectrum within the time resolution of 
150fs and its decay within about 10 psiii 

In this paper, we clarify how the photoinduced CT 
excitation occurs in the ID dimerized Mott insulator, 
by employing the ID dimerized Hubbard model cou- 
pled with a uniform lattice distortion. By the exact 
diagonalization method and a perturbative calculation, 
we demonstrate that an interdimer CT state, which is 
the lowest photoexcited state, destabilizes the dimerized 
phase. The interdimer CT state is found to be super- 
position of two dispersive elementally excitations. We 
investigate the origin of the midgap state observed ex- 
perimentally in K-TCNQ iil Candidates for the midgap 
state are shown to be intradimer and interdimer CT exci- 
tations from the photoexcited state (not from the ground 
state). Relevance to K-TCNQ is also discussed. 



II. MODEL 
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FIG. 1: Dimerized ID system with alternating transfer inte- 
grals. The thick lines correspond to the larger transfer inte- 
grals, and the thin lines to the smaller ones. 

To treat the dimerized phase of the ID strongly cor- 
related electron system coupled with the lattice system, 
we use a model whose Hamiltonian is given by 

H = H e + —P8 2 , (1) 
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where TL C is for the electronic part and the second term 
gives the energy of the lattice system. The elastic con- 
stant of the lattice is denoted by /3, and 8 (0 < 8 < 1) 
is the lattice dimerization parameter. The system size is 
denoted by N. Taking the adiabatic approximation for 
the lattice system, we omit the kinetic energy term of 
the lattice. The electronic part 7i e is the ID dimerized 
Hubbard model defined by 

N-l 

u e = - + (-^4+1,^1,, +4,^1+1,-) 

a 1=0 

+ U^ni A m^ (2) 
i 



(q. ct ) is the creation (annihilation) operator 



-i,a 



where cj 

of an electron with spin a at site I, ni t , 
The nearest-neighbor transfer integrals alternate as t(l + 
5), t(l — 5), ■ ■ ■ , as shown in Fig.^ The parameter t gives 
the transfer integral of the uniform lattice, and is used as 
the unit of energy in this paper. The filling of electrons 
is set to the half and the boundary condition is set to 
be periodic. The Coulomb interaction U/t of K-TCNQ 
has been estimated to be about 6 — !8,2ii2Si2& Almost all 
results shown in our paper are obtained for U /t = 15. 



III. NUMERICAL METHOD 

To search important photoexcited states, we calculate 
the optical conductivity spectrum of the ground state 
l^o) given by2I 



a(u>) = D8(u) +a ICS (uj), 
where D is its Drude weight defined by 

7T 



D 



N 



2tt^ K<A„Wo>| 2 

n>0 



N 



E n — Eq 



(3) 



(4) 



with K being the kinetic term of the Hamiltonian @ , J 
the current operator defined by 



j= it j2[i + (-1)^(4+!,^,, 



c L c *+i,^)> ( 5 ) 



Eq the ground state energy, \ip n ) the n-th excited state, 
and E n the corresponding energy. The lattice constant 
is set to be unity in this work. The regular component 
<7 rcg (w) is defined by 
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UJ + it + Eq - n 



(6) 



where e yields finite broadening and is set at O.li in our 
calculations below. We carry out the full diagonalization 
to obtain precise excited states for a{u>) and other quan- 
tities. Thus the system size N is restricted to a relatively 
small value of N — 8. 



IV. NUMERICAL RESULTS 



A. Optical conductivity 
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FIG. 2: Optical conductivity spectra in the ID dimerized 
Hubbard model for U/t = 10 and 15. The insets show results 
of nearly decoupled systems (5 — 0.95). 

Figure|21shows the optical conductivity spectra <t(lu) of 
the model @ for U = 10 and 15. In both cases, <r{ui) has 
similar characteristics mentioned below. At S = 0, there 
are a few states spread over U— At < u < U+At. It is well 
known that they consist of holon-doublon pair excita- 
tions, forming the so-called holon-doublon continuum in 
the N — > oo limiti 28 i 29 i 30 With finite 8 introduced, there 
appear several peaks with small spectral intensity, sug- 
gesting that new optically-allowed states are generated 
by breaking of the translation symmetry. Accordingly, 
the photoexcited states at 8 = get weakened spectral 
intensity and some of them eventually disappear. How- 
ever, the lowest photoexcited state survives up to the de- 
coupling limit (see the insets of Fig. - We focus on this 
state as a typical example of the holon-doublon states, 
and call this \i/ja) m the following. 

As 8 is further increased (8 > 0.4), a new optical state, 
which is called \i/jb) hereafter, appears as an exciton-likc 
sharp peak slightly above u> = U . The appearance of 
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this sharp peak is predicted in several preceding stud- 
iesi 26 i 31 i 32 i 33 i 34 This exciton-like state is found to belong 
to a different class from the holon-doublon continuum. 
Thus we focus on the two states, |^>a) and IV'b), as typ- 
ical photoexcited states in the dimerized system. The 
nature of these two states is discussed in the following. 



B. Adiabatic potentials 

In this subsection, we show the adiabatic potentials of 
the relevant states, \ipo), IV'a) and |Vb), to discuss their 
spin-Peierls stability. The adiabatic potential Ei(5) is 
given by 

Ei(S) = E?(S) + ^V, (7) 

where Ef(8) denotes the energy of the eigenstate \ipi) of 
H c {i=0, A, or B). Figure |3] shows the results for U — 
15. The adiabatic potential of is plotted only in 
the region where IV'b) has a pronounced spectral peak 
in (t(lu). The elastic coupling j3/t is set to 0.6 here for 
the adiabatic potential of |i/'o) to have a minimum at a 
considerably large value of the dimerization [6 ~ 0.6) to 
help us observe dimerization effects clearly. 

We have found that £a(<5) has a minimum at 8 = 0, 
which demonstrates that the dimerized state is destabi- 
lized by the optical excitation from to IV'a)- Hence, 
when the state IV'a) is photoexcited, the system can move 
toward the uniform state (6 = 0) along the adiabatic po- 
tential Ea(5) unless the recombination of the holon and 
the doublon occurs; this is the realization of the inverse 
spin-Peierls transition. By contrast, Eb(S) has a mini- 
mum around 8 — 0.4. This suggests that \i/jb) does not 
annihilate the lattice dimerization although it can weaken 
it. 



V. DECOUPLING LIMIT 

In this section, we present a simple picture of the rel- 
evant photoexcited states |Va) and |Vb) to understand 
their nature. For this purpose we introduce a perturba- 
tion theory that is justified in the decoupling limit 5 — ► 1. 
Since both of \ipA) and \?Pb} have finite spectrum inten- 
sity in the nearly decoupled system, it is reasonable to 
deal with both states from the decoupling limit, where 
the perturbative treatment of interdimer interactions is 
valid. 

Several studies have presented different approxima- 
tions for the ID dimerized Hubbard model to discuss the 
characteristics of ct(cj), for example, the dominant spec- 
tral intensity around u ~ £/,26 i 32 4 33 J 34 However, these 
preceding theories are justified only in the U/t — > oo 
limit. By contrast, our approximation can treat U/t ex- 
actly, and is proved to be appropriate to construct a quite 
simple picture of the photoexcited states. 




FIG. 3: Adiabatic potentials Ei(5) of three states (i =0, A, 
and B) for the N — 8 system. 



A. Hamiltonian 

In the decoupling limit (5 — > 1), the electronic part Ti e 
can be categorized into the intradimer part 7io and the 
interdimer part TLi as follows: 

H e = TCo + TCi, (8) 

with 

+U(ni,ini > i + ni+i,'[ni + i t i)] (9) 

and 

Hi = -t(l - S) £X)(4fl, < r<V + 4,aCl+h*), (10) 

where X)zec(o) re f ers to the summation over even (odd) I. 
Since |1 — 8\ « 1, Ho is the unperturbed Hamiltonian, 
and H\ is the perturbation term. The current operator 
is also rewritten as follows: 

J=J + Ji, (11) 

where 

Jo = it(l + S) Y,Y,( c t* c l+h° - (I 2 ) 

a Zee 

and 

■h = it(l - S) ^^(c, r /i +v - 4 +l a c LrT ). (13) 

a leo 
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B. Ground state 



Then its excitation energy is given by 



Let us first consider the Oth-order ground state of the 
system. In the completely decoupled case (5 = 1), the 
system is equivalent to Nd(= N/2) isolated dimers, each 
of which is described by the 2-site Hubbard model. Thus 
all the eigenstates are represented as the direct products 
of the eigenstates in the 2-site Hubbard model, whose 
Hamiltonian, eigenstates, and eigenvalues are summa- 
rized in Appendix. 

The Oth-order ground state |^q) is the direct product 
given by 



E 



|$) - |G) ® |G)i®---® |G) 



A ! ,/ 



(14) 



where \G) n denotes the ground state of the n-th dimer 
(n = 0, 1, • • • ,Nd~ 1). Note on the properties of \G) n in 
the U/t — > oo limit that, as shown in Appendix, \G) n is 
equivalent to the singlet dimer state of two S=l/2 quan- 
tum spins for U/t — > oo. Thus |V>o) m this limit is the 
product of the singlet dimers. In our case with finite U, 
IV'o) contains a small component of doubly occupied site. 
The Oth-order ground state energy is given by 



E° 



^e-[t(l + 5)\ 

N 
~2 



£7/2- y/U 2 /4 + 4t 2 (l + 5)' 



(15) 



where <T it) is the ground state energy of the 2-site Hub- 
bard model with a transfer integral t. It should be noted 
that the energy (|T3|) is exact to the first order in Hi 
because of 



(r \n 1 \r ) = o- 



(16) 



Intradimer CT state 



One class of relevant photoexcited states in the decou- 
pling limit consists of intradimer charge-transfer (CT) 
states, which are obtained by operating Jo onto IV'o)- 
We here introduce the following basis set: 



||0 2 ))„ = |G) ®---® K?>n-i® \0 2 ) 

® \G) n +x®--®\G) Nd - X) 



(17) 



where \0 2 ) n is the intradimer CT state of dimer n [see 
Eq. l|X27|) ]. Then we obtain 



N d -1 



JolVo ) =2it(l + 6)a ]T ||0 2 »„. 



(18) 



n=0 



The right-hand side of Eq. I|18fl is a linear combination 
of intradimer CT excitations IIO 2 ))^, and all the states 
are degenerate eigenstates of Ho with energy 



F° - I - - 1 

-^intra — I 2 



■[t(l + S)] + U. 



(19) 



u/2 + y/mj4 + W{TTW. 



(20) 



We note that these formulae for E® ntY& and AP ntra are also 
exact to the 1st order in Hi because of 



((0 2 \\ n Hi\\0 2 )) 7 



(21) 



for any n and m. Equation l|21|) implies that the local CT 
excitation \0 2 ) n transfers to the nearest-neighbor dimers 
(n ± 1) by at least 2-nd order processes of Hi. 



D. Interdimer CT excitations 

1. zeroth order in Tii 

The other class of important photoexcited states con- 
sists of interdimer CT excitations. Before discussing the 
details of these states, we introduce another useful basis. 
Let us consider a state where the ni-th dimer is in \X l a ) 
and the n2-th (rii < 712) dimer is in \Yg~ l ). Both \X l a ) 
and \Y^~ l ) are eigenstates of the 2-site Hubbard model 
and the superscripts, (1,4:— I) = (1,3) or (3,1), denote 
the number of electrons in each dimer. The spin variable 
a takes f or J., and a denotes the opposite direction of a. 
We assume that the other dimers are in the ground state 
\G). Such a state is represented by 



= |G)o®"-®|G>ni-l®|^i)ni 

® |G)„ 1+ i® •••® |G)„ 2 _i 

® \Yt l )n 2 ® \G) n2 + i ® • • • 

® \G) Nd -i. (22) 



The interdimer CT excitations are generated by oper- 
ating Ji onto l^o) • Several terms thus created are classi- 
fied into the lowest-energy terms and higher-energy terms 
as follows: 



Ji\^° ) = it(l-S) 



(a + /3) 



2 N d-i 



+ 1 



n=0 
l3\ 



\0\E^) n ^ n+ i — \E^O^/ n , n +i 
\E\0^) n ^n+l) 

higher energy terms. 



(23) 



All the states in the first part of Eq. I|23|) are degenerate 
eigenstates of Ho with energy 



F? =lfL 

inter \ 2 



[t(l + S)] + U-2t(l + S). (24) 



Its optical gap is given by 



^intcr — -^intcr ^0 



-2t(l + S) + 2^U 2 /4 + U 2 (l + S) 2 . (25) 
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FIG. 5: Intradimer and interdimer CT states, and their exci- 
tation energies for U/t >> 1. The spin indices in \Ot) and in 
\Ea) are omitted. 



FIG. 4: Schematic picture of the ground state, the intradimer 
and the interdimer CT states. The ellipses show the singlet 
dimer states, and the up (down) arrows are electrons with up 
(down) spins. 



It is obvious that AP ntra > A? ntcr for U > 0. Therefore 
the interdimer CT state is the lowest optical excitation. 
This conclusion can be easily understood in a large U 
system, which is appropriate to our study. For U/t >> 1, 
the excitation energies (1201) and l)25|) are evaluated as 



o 

intra 


inter 



u 
u 



4(1- 
2t(l 



6) 2 t 2 /U, 
\-S) + 8(1 



5ft 2 /U. 



(26) 
(27) 



The Coulomb term U appears because of the creation 
of double occupation by the charge transfer (see Fig. |3}. 
The last terms of Eqs. <|26[) and i|27|) originate from the 
destruction of the singlet dimer(s). The most important 
is the second term of Eq. J23): this is caused by the 
kinetic energy gain of the electron in \E*) and that of the 
doubly occupied site in \0%). That is, in an interdimer 
CT state \E^Og), the kinetic effect caused by moving 
the electron and the doubly occupied site reduces the 
excitation energy. 

Now, it becomes clear which CT state corresponds to 
IV'A(B))- Because the energy of |?/>a) is lower than that of 
|-0b) , IV'a) is the interdimer CT state and the |V>b) is the 
intradimer CT state (see Fig. This conclusion is quite 
natural because the spectral intensity of |$b) becomes 
dominant while that of \ipA) decreases as S approaches 1. 
In addition, our interpretation of |^b) is consistent with 
that in the preceding studiesi 26 i 31 i 32 i 33 i 34 



the lst-order degenerate perturbation theory to obtain 
the energy gap to the 1st order in Hi. For this purpose 
we calculate the matrix elements of Hi for the degener- 
ate states \X l a Y^ l ) ni , n2 , with (X l ,Y 4 ~ l ) = {0 3 ,E r ) or 
(E 1 ^ 3 ), to obtain 



Hi\X l a Y 4 ) ni ,„ 2 



t(l-S) 



{a + p) 2 (\X l a Yt l )n lM 



X a Y„ ) ni , 7i2-i + \X l a Yg ) ni +i, na 

X a Y s )m-l,ri 2 ) 

for n 2 - ni > 1, (28) 



and 



HAXlYf 



t(l-S) 



(a + P) 2 {\X l a Yt l ) ni ^ + 



a^a /m,«2 



+ \ X l Y & )m-l,n 2 ) 
for n2 — n± = X. 



(29) 



Here, other states with higher (or lower) energies are 
omitted because they are not needed at present. 
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2. First order in Hi 



Here we discuss the lst-order correction to the opti- 
cal gap due to the interdimer CT state. We have found 
that the Oth-order interdimer CT state l|23[) consists of 
a number of interdimer CT excitations \0 3 E^} n ^ n+ i and 
|^0|)„,„+i (0 < n < N d ). All these CT states have 
degenerate energy l)24[l. In addition, states that contain 
\0 3 ) and \E\) in separate dimers rti,rt 2 i\n% — ri\\ > 1) 
also have the same energy 1)240 ■ Thus we need to perform 



FIG. 6: Dynamics of excited dimers in the optically excited 
states. Effective transfer integrals for |O s ) and l-E 1 ), and that 
for \0 2 ) are very different as indicated. 

The relations l|28ll and l|29|l show that the local excita- 
tions \0 3 ) and \Eg) move in the system by the lst-order 
process of Hi, as illustrated in Fig. EI This is quite dif- 
ferent from the intradimer CT state, where the local CT 
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excitation \0 2 ) transfers by at least 2nd-order processes 
ofW x . 

The lst-order Hamiltonian Hi defined by and Hty 
can be diagonalized analytically in the subspace spanned 
by \X l a Yg~ l ) nitn2 . We here introduce a new basis set: 



f 



4°l)™i>™2 - \EiO^) ni>n2 ] . (30) 
By using this basis set, Ji\iPq) is written as 

2 ~ 



|nx,n 2 ) = - L|Of^)„ 1>n2 + lOJ^)^^ 



Jrl^o ) = it{l-S)^-BL £ \n,n + l), (31) 



n=0 



where the higher-energy terms are omitted. Then we 
obtain 

Hi\n u n 2 ) = ^ 7 ^ (a + /?) 2 (K + 1, n 2 ) + |m - 1, n 2 ) 



+|ni,n 2 + 1) + |ni,n 2 - 1)) 
for n 2 — ni > 2, 



(32) 



and 



WiK,n 2 ) = ^ - S \ a + f3) 2 (\ ni - l,n 2 ) + K,n 2 + 1)) 



for n 2 — ni = 1. 



(33) 



It is noted that the following boundary condition should 
be used: 



|m,n 2 = N d ) = |0,m) 

|m = -l,n 2 ) = \n 2 ,N d -l). 



(34) 



Then this problem is reduced to the spinless free-fermion 
system with two particles under the antiperiodic bound- 
ary condition. The eigenstates are written as 



Mfci,Aa)) - 1/N d l n i' n2 ) 

0<rii<n 2 <A r d-l 



X e 



(35) 



where 



«1(2) 



2m 1(2) + 1 



7T (m 1(2) =0,1,-- - ,AT d -l), (36) 



with fci ^ /c 2 . Their energies are given by 



E\k u k 2 



-(a + /3)'[cos(fc 1 )+cos(fc 2 )]. (37) 



Since \ip%) is the lowest-energy state of Eq. (|35f) with 
the total wave number being zero, this state is given by 

N — 2 

h&) = M*,-*)) with k = ^-ir. (38) 



For jV = 8, the wave numbers are ki = 3/47T, fc 2 
— 3/47T, and then we have 

El ntcl .(k 1 ,k 2 ) = - t -^l(a + pr. 

Thus the optical gap to the 1st order is 



(39) 



A 1 _ pO 

inter — inter 



^inter 



El 



(40) 



In Fig. El the perturbative results for the optical ex- 
citation energies A? ntra , AP ntcr , and A? ntra are compared 
with the exact results obtained by the numerical diago- 
nalization. The numerically calculated, optical excitation 
energies of IV'A(B)) a re denoted by Aa(b)- The analytical 
results for AP ntra show good agreement with the numeri- 
cal results for Ab- Hence we confirm that the state |^b) 
is the intradimer CT state. The results for A- ) ntcr agree 
with the numerical results for Aa around 5 ~ 1. As S de- 
creases, the difference between them becomes relatively 
large compared to that between A? ntra and A B because 
of the finite lst-order correction due to the dynamics of 
\0^) and \Eg). The agreement between A i 1 ntor and the 
diagonalization results supports our argument on the lst- 
order perturbation. 




08 09" 



FIG. 7: Optical excitation energies. The symbols show the 
exact-diagonalization results, and the lines the perturbative 
results. 



3. Origin of the "inverse" spin-Peierls instability 

Up to now, we have clarified how the CT occurs in 
the relevant photoexcited states in the dimerized phase. 
Here, we focus on the relation between the CT type and 
the (in)stability against the lattice dimerization. The 
dimerization in the spin-Peierls system is caused by the 
energy gain due to the spin singlet states formed on the 
dimers. Hence whether the lattice dimerization is stable 
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or not in the CT states depends on how the spin singlet 
states are destroyed by the CT excitations. To examine 
this, we calculated a quantity that measures the strength 
of the spin dimcrization: 



with 



D = \(S ■ Si) - (5i • S 2 ) 



(41) 



The obtained results are shown in Fig. |H1 The ground 
state has large spin dimerization even when S is not so 
large: at 6 — 0.2, D is about 80% of its maximum value. 
By contrast, that of |t/a) is strongly reduced, which in- 
dicates that a large part of the spin dimers is destroyed 
by the interdimer CT excitation. In IV'b), D has a larger 
value than that of \4>a) although it is of course smaller 
than that of \G). This result is consistent with the fact 
that the strength of \i/>a) f° r destabilizing the dimer phase 
is much larger than that of | V'b)- 
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ijj + ie + Ea - H 



(43) 



Nuj 



-Im 
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i 



uj + ie- E A +U 
Here, D x is the Drude weight of |Va) defined by 



2 TV 

N 



n^A 



E n — Ea 



(44) 



(45) 



where Ea is the energy of IV'a) 
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FIG. 8: Spin dimerization in the three states. 



VI. MIDGAP STATES IN ai(u) 
A. Numerical results 

Here, we discuss optical responses of the system after 
the photoexcitation. Our interest in this aspect has been 
stimulated by the experiment by Okamoto et ai, who 
have reported the appearance of a midgap state in the 
reflectivity spectrum of K-TCNQ immediately after the 
photoexcitation. 3 '' 1 ' Because the midgap state is observed 
after the photoexcitation causing the inverse spin-Peierls 
transition, the realized state is assumed to be IV'a)- In 
addition, we note that it is unnecessary to take account 
of the lattice relaxation because the midgap state already 
appears before the lattice starts to moveiii Then the 
key quantity is the optical conductivity in the state |^a) 
given by 



(42) 



FIG. 9: Optical conductivity spectra in the photoexcited state 
|^a) for U/t = 15. 

Figure [5] shows numerical diagonalization results of 
ctiH for U/t = 15 and N = 8. At 5 = 0, there 
appear the large Drude peak and another pronounced 
low-energy peak. The presence of the Drude peak sug- 
gests the metallic property caused by the photodoped 
carriers, the holon and the doubloniS 5 ^^ 6 . The latter low- 
energy peak corresponds to the excitation from to 
the nearly degenerate excited state (named | V'A' ) ) j25i2& 
which is an important feature of ID uniform Mott insu- 
lators^ As 5 increases, these peaks lose their intensity 
and several new peaks appear. The high-energy peak 
around to ~ U corresponds to the intradimer CT exci- 
tation from \G) to |0 2 ), which is also observed in <t(uj). 
Since I^a) still includes a number of dimers in the state 
\G) except for the charge-transfered dimers \E^O%.}, this 
peak appears. More interesting are midgap peaks around 
u> ~ At. The midgap peaks develop with increasing 5, 
indicating that they are closely related to the dimeriza- 
tion. At first glance, the midgap state seems to belong 
to a different class from the nearly degenerate excited 
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state IV'A')- The midgap state appears only for finite <5 
while \ipA_r) exists even for 5 = 0. In addition, the energy 
of the midgap state (oj/t ~ 4) is far from that of \il>A') 
(iv/t ~ 1). Since IV'A') is close to \iPa) even for finite 6, 
\ipA') would be a linear combination of |£^0?), which are 
also the constituents of \ipA)- By contrast, the midgap 
state consists of other types of excitations, as discussed 
in the next subsection. 



B. Zeroth order in Hi 

As a starting point to discuss the origin of the midgap 
states, let us present again the analytical form of I^a): 
the perturbative result of the wave function of \ipA) is 
given by Eq. Ij35(l with k\ = — k 2 = k, 



\o = VN d 



r ik(n 1 ~n 2 ) _ g — ik(m— n 2 )l 



n\ <ri2 



— \E^O^) niyn2 — |£/JOj) nii „ 2 ), 



(46) 



with k = ir(N — 2)/N. Since the spectral intensity of 
the midgap states increases with 5, candidates for the 
midgap states would be the intradimer CT excitations 
obtained by operating Jq onto \i[> A ) as follows: 



it(l + S) 
N d 


E 

111 <"2 


e ifc(ni -Tl2 


>-e 


x \(\EM 


ni ,n 2 


+ \eM) 


ni ,n 2 


+ \opl) ni 


,112 "T" 


o\o}) ni 


ri2 


- \E^Ef) ni 


,™2 _ 




n 2 


- \o]o\) ni 


,"2 _ 


oJo T 3 ) ni 


12)1 



-ik(nx— n 2 )l 



(47) 



where other states with higher energy are omitted. 
All the states in the above equation, \E l a E^ '•*) and 
|0^0^ 4 ~^), are degenerate eigenstates of Ho with energy, 



N 



E° midx = — [U/2 ^U 2 /A + At 2 (l + 5) 2 
+2 v /C/ 2 /4 + 4t 2 (l+<5) 2 . 



(48) 



In addition, we find that another class of states con- 
tribute to the midgap states. Operating Ji onto \ip A ), we 
obtain 

Ji\ti) =-^fc^(a 2 -/3 2 )sin(fc)^||i? 2 ))„, (49) 



where the new state \ \E 2 )) n is defined by 

\\E 2 )) n = \G) ®---®\G) n . 1 ®\E 2 
<g> |G) n+ i ® • • • <g> \G) Nd -i. 



||i? 2 ))„ is also an degenerate eigenstate of Hq with en- 
ergy, 



£L, mid2 



N 



U/2- y/U 2 /4 + U 2 (l+5)' 



+2 v /[/ 2 /4 + 4< 2 (l + ( 5) 5 



— -^midl- 



(51) 
(52) 



Figure ITU1 displays the relevant states and their relative 
energies given by 



-^midl ^inter 
- C/ mid2 -^intra 



2t(l + 5), (53) 

e+[t(l + «)]-!/ 

4(1 + rJ) 2 t 2 /^ 

for U/t»l. (54) 



Here it should be noted that the energy difference l|53|) 
between |^) and the excited state JoIV'a) 1S quite close 
to the gap (~ At) observed in Fig. 03 Thus these excited 
states are considered to be the midgap states. 



4(i+5) 2 f 2 /U 



-|E 2 > 



|G> 



2t(l+8) 




= |o 2 > 




2 t 2 /U 


u- 




+ 



|0 3 >|0 1 > 
or 

|E3>| E 1> 

|Va> = |o3>|eO 

21(1+5) 
8(1+5? t 2 /U 

|G>|G> 



FIG. 10: Midgap and other relevant states. The energy dif- 
ferences are also shown for U/t » 1. 



C. First order in Hi 

Without perturbation, the midgap states are degener- 
ate. However, the numerical results in Fig. [!|] show that 
these states are split for finite (1 — 6) by the perturba- 
tion. Then we discuss the effect of the interdimer in- 
teraction Hi on the midgap states to the 1st order. The 
following degenerate states should be taken into account: 
\ElEl) nun2 , \OlOl) nun2 , and ||£: 2 )) n . For simplicity, 
we introduce a basis set given by 

|nI7W) = \(\E^E\) nun2 + \ElE]) ni ^ 

- |Ot°?>*i.*« ~ l i°T>«i>»2) for "1 < "2 



(50) 



■^(\E^Ef) n2tni + \E\E^) n2 ^ ni 
"tODna.m - \0\0)) m ) for n 2 



2 

10 



< m. 

(55) 
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By using this, we represent the relation l|47|) in a simpler 
form: 



MO = it{1 + s) 



N d 

ni<U2 

x (|ni,n 2 ) - \n 2 ,ni)) 



y ' Jgifc("i —n 2 ) _ g— *fc(ni— n 2 )j 

(56) 



The matrix elements of 7ix are calculated as follows: 
t(l 8) 



Hi\n 1 ,n 2 ) = 



4 

t(l-6) 



{a + /3) 2 (\m,n 2 + 1) + \m,n 2 - 1)) 

(a - /3) 2 (|ni + l,n 2 ) + K - l,n 2 )) 

(57) 



for n 2 — ni > 1, 



HilnTTW) = ^— -^(a + /3) 2 |m,n 2 + 1) 



4 

*(!-*) 
4 

*(!-*) 
2 

t(l-5) 



(a-/3) 2 |m-l,n 2 )) 
{a + (3) 2 \\E 2 )) ni 



(a-/3) 2 \\E 2 )) n2 (58) 



for n 2 = ni + 1, 



Hilm.na) = — ; ^ (a + /?) 2 |ni,n 2 - 1) 



4 

*(!-*) 
4 

*(!-*) 
2 

t(l - 5) 



(a-/?) 2 |m + l,n a » 
{a + (3f\\E 2 )) ni 



(a-(3) 2 \\E 2 )) n2 (59) 



for n 2 = iii — 1) and 

t(l-*) 



(a + /?) 2 (|n,n + l)- |rc,ra-l)) 



!i^_i) (a _ / 3 ) 2 ( | 71 + 1;n) _| n _ ljn)) . 



(60) 

Equations l|57|) , 1581) , (|59l) , and i|60[l show several charac- 
teristics of the local excitations included in the midgap 
states. For example, the local states \E^) and \OP) 
move in the system with the larger transfer integral 
-t(l - <5)(a + /3) 2 /4 while the states \O x a ) and \Ef) move 
with the smaller transfer integral t(l — 8) (a — /3) 2 /4 (see 
Fig. Illfl . In particular, if the terms containing ||i? 2 )) 
are omitted, the system becomes equivalent to the spin- 
less free-fermion system with a light "particle" (\E^) or 
\Ol)) and a heavy "particle" (|C£) or The terms 

containing ||-E 2 )) mean that a collision between the light 
"particle" and the heavy "particle" creates ||i? 2 ))- It is 
also found that the state Hi? 2 )) splits into a pair of the 



light "particle" and the heavy "particle". Here we note 
that the following boundary conditions should be used: 

\ni,n 2 = N d ) = -|ni,0) 

\ni,n 2 = -I) = -\ni,N d -l) 

| Til = N d ,n 2 ) = -|0,n 2 ) 

|ni = -l,n 2 ) = -\N d -l,n 2 ) 

\\E 2 )) Nd = \\E 2 )) 



\E 2 ))-i 



\E 2 )) Nd - 



(61) 



t(l-S)(a-p) 2 /4 
-|E 3 > 



JG> m m |G> 



| E 1>. 



-t (l-8)(a+pf/4 



FIG. 11: Dynamics of excited dimers in the midgap state. It 
contains a light "particle" IE 1 } and a heavy "particle" \E 3 }. 

The effective Hamiltonian derived here is difficult to 
diagonalize analytically. Hence we carry out the numer- 
ical diagonalization to obtain the i-th eigenstate in the 
form: 



\ti 



+ 



E 

0<ni<n2<A r d-l 

E 

0<ri 2 <ni<A r ti -l 

E 

0<n<7V rf -l 



a l {ni,n 2 )\ni,n 2 ) 
a l (ni,n 2 )|ni,n 2 ) 



6'(n)||£; 2 )) r 



(62) 



To pick up optically allowed excited states, we calculate 
the spectral intensity to the 0th order in (1 —8): 



(63) 



where 



/,o,Ti/i \ 2t(l + <J) 

m\JoWmid) = Nd 



x {E sm[k(n 2 - m)]a l (ni,n 2 ) 
+ ^2 sin[k(n 2 - ni)]a i (ni,n 2 )}(64) 

n 2 <n 1 

In Fig. the perturbative results for the energies 
and the spectral intensities are compared with the exact- 
diagonalization results for N = 8. The numerical results 
show that there are three photoexcited states, which are 
well reproduced by our perturbation theory. We obtain 
two isolated states, IV'iid) and IV'mid)' an< ^ several de- 
generate states that are collectively named IV'mid)- The 
spectral intensity of | - 0^ lid ) is obtained by the summation 
of those of the degenerate states. The higher-order effect 
would lift the degeneracy for IV'mid)- Some of thus split 
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states must be optically active. As for the spectral inten- 
sity, the perturbative results are good only in the vicinity 
of 5 = 1. This is because the obtained wave function is 
correct only to the 0-th order in 1 — 5 and also because 
the interdimer current operator J\ is omitted in Eq. (|63ll . 



14 



CD 



13.5 



13 



(a) 



U/t=15 

^\ on Exact 

A diagonalizatiorj 



X 



V 3 

\ □ v / 

□ XlV mid> 



0.6 




8 0.8 



c 

CD 



22 
o 

CD 
CL 
00 







(bj 1st-order 

perturbation 




8 0.8 



FIG. 12: Calculated energies and spectral intensities for the 
midgap states. The symbols show the results of the exact 
diagonalization of the Hamiltonian (|5J and the lines denote 
the results of the lst-order perturbation theory. The circles 
and the solid lines are for I'i/'mid)! the squares and the dotted 
lines for [V'mid)) an d the triangles and the dashed lines for 

h/4id>- 



VII. SUMMARY AND DISCUSSION 

We have studied the properties of the photoexcited 
states in the ID dimerized Hubbard model coupled with 
the lattice. A notable point of our conclusion is that the 
lowest photoexcited state in this model is the interdimer 
CT state. It is true that this point is inconsistent with 
the experimental result for K-TCNQ, where the lowest- 
energy peak in a(uj) of K-TCNQ has dominant spectral 
intensityiSiiSiiS This experimental fact implies that the 
lowest photoexcited state is an intradimer CT state, as 
the authors pointed outi^&iS This discrepancy would be 
due to the finite nearest-neighbor Coulomb interaction, 
often denoted by V, of K-TCNQ. The intradimer CT 



state could lower its energy by the effective attractive 
interaction V between the holon and the doublon. This 
is because, in the intradimer CT state, the mean distance 
between them is smaller than that in the interdimer CT 
state. The expected energy diagram for a finite- V system 
is compared with our results for V — in Fig. 1131 

However, in the experiments of the photoinduced in- 
verse spin-Peierls transition, the excitation energy so far 
corresponds to the energy of the interdimer CT state 
rather than that of the lowest-energy intradimer CT 
state. 10 ' 11 In this sense, the experimental result is con- 
sistent with our conclusion that the interdimer CT state 
destabilizes the dimerized phase, leading to the inverse 
spin-Peierls transition. Our assumption that the realized 
state after the photoexcitation is the interdimer CT state 
is justified in discussing the photoinduced dynamics ob- 
served so far. If an experiment were performed in which 
the energy of the pulsed laser is lowered to coincide with 
the excitation energy of the lowest-energy intradimer CT 
state, the dimerization would be only weakly reduced 
and the following lattice dynamics would reflect the lo- 
cal character of the intradimer CT state, which is quite 
different from the character of the interdimer CT states. 

The photoconductivity experiment of K-TCNQ also 
supports our conclusion on the photoexcited state. The 
experimental results shown in Fig. 2 (a) of Ref . ITTI reveal 
the conducting nature of the photoexcited state. With 
the free-carrier-like charged dimer states l-E 1 ) and |0 3 ), 
the interdimer CT state is the best candidate for the pho- 
toexcited state. By contrast the intradimer CT state has 
an exciton-like bound state that does not have the con- 
ducting nature. This situation would remain unchanged 
for finite V systems. 

Here we will give a comment on another way to de- 
scribe the photoexcited states of the ID dimerized Mott 
insulators. It is the odd (or even) CT state proposed in 
Ref. [3^. Since the dimerization occurs, a CT state with 
even parity at 5 = becomes optically active acquiring 
a small spectral weight^ This description is valid for 
small-5 systems. It is a complementary approach to our 
treatment. For example, the interdimer CT state adia- 
batically connects to the lowest photoexcited state of the 
uniform system, which corresponds to the odd CT state. 
Both the interdimer and intradimer CT states are linear 
combinations of the odd and even CT states. As S in- 
creases, the degree of hybridization changes continuously. 
Hence, our discussion from the decoupling limit, which 
allows us a relatively simple analytical treatment, is help- 
ful to understand the optical excitations of the present 
system. 

Another conclusion of our study is on the origin of 
the midgap state observed just after the photoexcitation. 
Our calculation has clarified that the midgap states cor- 
respond to both the intradimer and interdimer CT ex- 
citations from \ipA}- Here we note that the degeneracy 
between \ \E 2 )) and ^O 1 ) or l^ 3 ^ 1 ) would be lifted for 
finite V (see Fig. 11311 . Then the peaks might be split 
into two parts. In any case, a dominant spectral weight 
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belongs to the intradimer CT excitation due to the large 
transfer integral, not to the interdimer CT excitation. 
Therefore the observed midgap state in K-TCNQ would 
be due to the intradimer CT excitation from the pho- 
toexcited state IV'a)- 

Here we discuss another possibility that the experimen- 
tally observed midgap state is not caused by the purely 
electronic mechanism. Then the most probable candi- 
date for the midgap state would be a polaron where a 
photodoped carrier is coupled with phonons via electron- 
lattice couplings. In fact, Okamoto et al. have found 
that the photoinduccd phenomena in K-TCNQ are af- 
fected by lattice oscillations due to three relevant phonon 
modes. 11 However, it seems inappropriate to explain the 
midgap state by those phonon modes. The lowest en- 
ergy (20cm _1 ) phonon mode has a too- long oscillation 
period, and the other higher-energy modes decay much 
faster than the midgap state (which decays in about 3ps). 
In this paper, we have provided an idea only on the gener- 
ation of the midgap state, taking account of its ultrafast 
development. It is a challenging problem to discuss why 
the midgap state has a relatively long decay time (3ps), 
although it is beyond the scope of the present study. A 
possible scenario is that the midgap state caused by the 
purely electronic mechanism couples with slowly decay- 
ing phonon modes, such as the 20cm~ 1 mode. 



(a) V=0 (our results) 



E 3 >|E 1 > 




= |0 3 >|E 1 > 
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APPENDIX: 2-SITE HUBBARD MODEL 



The 2-site Hubbard model has the Hamiltonian given 



h= -i^(c| CT c 0(T + c ( ! )(T ci (T ) + [/(no T nox+niT«ii), (A.l) 



and the current operator defined by 

3 = itj^^o^ - c L c 0o 



(A.2) 



Below we summarize the eigenstates, their energies, and 
the matrix elements of the current operator between 
them, for different numbers of electrons. 



1. 1-electron states 

When the dimer contains only one electron, the eigen- 
states of the Hamiltonian i|A.l|) for t > are: 

\El) = -L( c t CT + 4)|0) (P = l), (A.3) 

Pi) = -±=(4-4)10) (P = -l). (A.4) 

In our notation, E(0) indicates that the state has the 
parity P = 1(— 1), the superscript shows the number 
of electrons, and the subscript a is the spin index for 
the dimer or an electron. These states obey the eigen 
equations: 



h\El) = -t\El), 
h\Ol) = t\Ol). 



(A.5) 
(A.6) 



FIG. 13: Schematic picture of the relevant optical excitations 
for (a) V = and (b) finite V. 



The optical properties are given by 

j\El) = it\Ol), 
j\Ol) = -it\El). 



(A.7) 
(A.8) 
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3-electron states 



where 



Next, we consider the case with three electrons. The 
eigenstates are 

\El) = -1(444-444)10), (A.9) 

\od = i=(4cL4 + 444)10), (a.io) 

where a =\ (J.) for a =\ (|). The eigen equations and 
the optical properties are summarized as 



h\El) = (U + t)\E 3 a ), 
h\Ol) = {U-t)\Ol), 



and 



j\El) 



-it\Ot 



j\Ol) = it\El). 



(A.ll) 
(A.12) 



(A.13) 
(A.14) 



3. 2-electron spin-singlet states 

Finally, we treat the states containing one electron 
with up spin and one electron with down spin. In this 
subspace, relevant states to photoexcitation are singlet 
eigenstates summarized as 



\o 2 ) = 

\G) = 
\E 2 ) = 



V2 

a\d)+(3\c), 
-j3\d) + a\c) 



(A.15) 

(A.16) 
(A.17) 



and 



Q 
P = 

c = 



\c) 
\d) = 

= 2t/C, 



— (44 + c-^CqjJIO), 



y/lP/± + U 2 -U/2 jC, 



U 2 /2 + 8t 2 - Uy/U 2 /^ + U 2 



1/2 



(A.18) 
(A.19) 



(A.20) 
(A.21) 

(A.22) 



We have a — *■ 1 and (3 — > in the U/t — ► oo limit. Then 
the ground state \G) approaches \d), which is equivalent 
to the so-called dimer singlet. Eigen equations are sum- 
marized as follows: 



h\0 2 ) = Up 2 ), 
h\G) = e-(t)\G), 
h\E 2 ) = e+(t)\E 2 ). 



where 



~-(t) = U/2±y/U 2 /4 + 4t 2 . 



The optical property is given by 



j\G) = 2ita\0 2 ), 
j\0 2 ) = -2V2it\d) 

= -2V2it{a\G) ~P\E 2 )). 



(A.23) 
(A.24) 
(A.25) 

(A.26) 

(A.27) 
(A.28) 



* Electronic address: maeshima@ims.ac.jp 

URL: http : //magellan. ims . ac . jp/~maeshima/ 

1 G. Yu, C. H. Lee, A. J. Heeger, N. Herron, and E. M. 
McCarron, Phys. Rev. Lett. 67, 2581 (1991). 

2 M. Fiebig, K. Miyano, Y. Tomioka, and Y. Tokura, Science 
280, 1925 (1998). 

3 S. Iwai, M. Ono, A. Maeda, H. Matsuzaki, H. Kishida, 
H. Okamoto, and Y. Tokura, Phys. Rev. Lett. 91, 057401 
(2003). 

4 S. Koshihara, Y. Tokura, T. Mitani, G. Saito, and T. Koda, 
Phys. Rev. B 42, 6853 (1990). 

5 S. Koshihara, Y. Tokura, K. Takede, and T. Koda, Phys. 
Rev. Lett. 68, 1148 (1992). 

6 H. Matsuzaki, T. Matsuoka, H. Kishida, K. Takizawa, H. 
Miyasaka, K. Sugiura, M. Yamashita, and H. Okamoto, 
Phys. Rev. Lett. 90, 046401 (2003). 

7 H. Matsuzaki, W. Fujita, K. Awaga, and H. Okamoto, 
Phys. Rev. Lett. 91, 017403 (2003). 

8 M. Chollet, L. Guerin, N. Uchida, S. Fukaya, H. Shi- 
moda, T. Ishikawa, K. Matsuda, T. Hasegawa, A. Ota, 
H. Yamochi, G. Saito, R. Tazaki, S. Adachi, and S. Koshi- 
hara, Science 307, 86 (2005). 

9 N. Tajima, J. Fujisawa, N. Naka, T. Ishihara, R. Kato, Y. 



Nishio, and K. Kajita, J. Phys. Soc. Jpn. 74, 511 (2005). 

S. Koshihara, Y. Tokura, Y. Iwasa, and T. Koda, Phys. 
Rev. B 44, 431 (1991). 

1 H. Okamoto, K. Ikegami, T. Wakabayashi, Y. Ishige, J. 
Togo, H. Kishida, and H. Matsuzaki, Phys. Rev. Lett. 96, 
(2006) 037405. 

2 K. Nasu, Rep. Prog. Phys. 67, 1607 (2004). 

3 K. Yonemitsu and K. Nasu, J. Phys. Soc. Jpn. 75, 011008 
(2006). 

4 E. Hanamura and N. Nagaosa, J. Phys. Soc. Jpn. 56, 2080 
(1987). 

5 N. Nagaosa and T. Ogawa, Phys. Rev. B 39, 4472 (1989). 

6 K. Koshino and T. Ogawa, J. Phys. Soc. Jpn. 67, 2174 
(1998). 

7 P. Huai, H. Zheng, and K. Nasu, J. Phys. Soc. Jpn. 69, 
1788 (2000). 

8 K. Yonemitsu and N. Miyashita, Phys. Rev. B 68, 075113 
(2003). 

9 K. Iwano, Phys. Rev. B 70, 241102 (2004). 

K. Yonemitsu, Phys. Rev. B 73, 155120 (2006). 

1 J. G. Vegter, T. Hibma, and J. Kommandeur, Chem. Phys. 
Lett. 3, 427 (1969). 

2 N. Sakai, I. Shirotani, and S. Minomura, Bull. Chem. Soc. 



13 



Jpn. 45, 3321 (1972). 

H. Terauchi, Phys. Rev. B 17, 2446 (1978). 

K. Yakushi, T. Kusaka, and H. Kuroda, Chem. Phys. Lett. 

68, 139 (1979). 

M. Meneghetti, Phys. Rev. B 44, 8554 (1991). 

F. B. Gallagher and S. Mazumdar, Phys. Rev. B 56, 15025 

(1997). 

The formula for the optical conductivity is given, for ex- 
ample, by B. S. Shastry and B. Sutherland, Phys. Rev. 
Lett. 65, 243 (1990). 

R. M. Fye, M. J. Martins, D. J. Scalapino, J. Wagner, and 
W. Hanke, Phys. Rev. B 44, 6909(1991). 
W. Stephan and K. Penc, Phys. Rev. B 54, 17269 (1996). 
E. Jeckelmann, F. Gebhard, and F. H. L. Essler, Phys. 
Rev. Lett. 85, 3910 (2000). 

E. Y. Loh and D. K. Campbell, Synth. Met. 27, A499 
(1988). 



S. K. Lyo, Phys. Rev. B 18, 5835 (1978). 

F. Gebhard, K. Bott, M. Scheidler, P. Thomas, and S. W. 

Koch, Philos. Mag. B 75, 13 (1997). 

F. Gebhard, K. Bott, M. Scheidler, P. Thomas, and S. W. 
Koch, Philos. Mag. B 75, 47 (1997). 

N. Maeshima and K. Yonemitsu, J. Phys. Soc. Jpn. 74, 
2671 (2005). 

N. Maeshima and K. Yonemitsu, J. Phys: Conf. Ser. 21, 
183 (2005). 

Y. Mizuno, K. Tsutsui, T. Tohyama, and S. Maekawa, 
Phys. Rev. B 62, 4769 (2000). 

H. Okamoto, Y. Tokura, and T. Koda, Phys. Rev. B 36, 
3858 (1987). 

S. Yamaguchi, Y. Moritomo, and Y. Tokura, Phys. Rev. B 
48, 6654 (1993). 



